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Abstract 

We study the gravitational backreaction of the non-abelian gauge field on the gravity dual to a 2+1 p-wave 
superconductor. We observe that as in the p + ip system a second order phase transition exists between a 
superconducting and a normal state. Moreover, we conclude that, below the phase transition temperature 
T c the lowest free energy is achieved by the p-wave solution. In order to probe the solution, we compute the 
holographic entanglement entropy. For both p and p + ip systems the entanglement entropy satisfies an area 
law. For any given entangling surface, the p-wave superconductor has lower entanglement entropy. 



1 rarias@fisica.unlp.edu.ar 
2 peznacho@gmail.com 



1 Introduction 



The AdS/CFT correspondence [TJ [2] [3] in its original form relates a conformal field theory in d dimensions 
with type II string theory on AdSd+i- The power of the correspondence arises from the fact that it is a 
weak/strong coupling duality i.e. it relates the strong coupling regime of the field theory with the weak 
coupling regime of the string theory and viceversa. After the works [4j [5], the gauge/gravity conjecture 
begun to be an useful tool to study condensed matter physics. In particular, it has been applied to study 
strongly correlated condensed matter systems through the analysis of a semiclassical dual gravity theory (see 
[8J for a review) . 

In the present work we analyze the backreaction of the gravity dual to a p-wave superconductoj^l in 
3+1 dimensions [6] (see [9l [10] for a similar treatment in 4+1 dimensions). Along the way we rederive the 
backreaction of the colorful p + ip superconductors previously studied in [TT] . We use the prescription given 
in [T2"] [T3"l [14] to compute the entanglement entropy from the holographic point of view for both gravity 
duals. Similar computations of entanglement entropy in backgrounds duals to condensed matter systems 
can be found in [T5] |T6l |T7] [15] . 

p-wave superconductivity is a phase of matter produced when electrons with relative angular momentum 
j = 1 form Copper pairs and condense. In other words, the operator that condense is a vector, charged 
under a U(l) symmetry. This kind of superconductivity is supposed to originate from "strongly correlated" 
electrons and therefore the BCS theory is not the correct approach to study its microscopic dynamics. This 
phenomena is a challenge for theoretical physics, and due to the fundamental property of the gauge/gravity 
duality mentioned above one could envisage the study of such systems through their weak gravity dual. We 
are going to introduce the minimal ingredients that one needs on the gravity side in order to reproduce 
the dynamics of the superconductor, this kind of approach aims to reproduce the properties of a condensed 
matter system without trying to explain their microscopical origin. 

String theory embeddings of p-wave superconductors were studied in \T§\ [5D] [5TJ [35] [231 HI] ■ The necessary 
minimal ingredients on the bulk to have finite temperature, chemical potential and spontaneous symmetry 
breaking (SSB) are: a black hole geometry and a non-Abelian gauge field (25] [26]. The solutions we will 
consider are asymptotically AdS backgrounds with a SU(2) non-Abelian gauge field. The SSB is realized 
on the bulk side as a non-trivial asymptotics (hair) for the gauge field. The chemical potential and the SSB 
arise by turning on two independent directions inside the non-Abelian gauge group. The symmetry breaking 
occurs on the gravity side through the formation of a condensate outsides the horizon. 

The entanglement entropy (EE) between a subsystem A and it's complement B is the von Neumann 
entropy 

Sa = -Tr A {p A \np A ). (1) 

Here p A = Trig (p) is the density matrix obtained by tracing the density matrix of the whole system p over 
the B subsystem degrees of freedom. Roughly speaking S A measures how much information is hidden inside 

3 The tiny difference between a superconductor and a superfiuid arises in the fact that although both effects are produced 
by a spontaneously symmetry breaking, in the first case there is a local symmetry which is spontaneously broken while in the 
second case is a global symmetry. We are going to use the terms superfiuid and superconductor interchangably here. For the 
considered phenomena this distinction does not make any difference. 
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B when we subdivide the system. From the point of view of the dual gravity theory the EE was conjectured 
\\2\ to be proportional to the bulk minimal area surface, 7.4, whose boundary at infinity coincides with the 
boundary of A (see [2] for a review) 

Sa = - 2 ■ (2 

Here k is the bulk gravitational constant. Note that the standard thermal entropy is obtained as a particular 
case of the EE, when the region A is the whole system. In [27J the authors provide a demonstration of 
this holographic technique to compute the entanglement entropy for spherical surfaces and zero temperature 
CFTs. On this work we compute this quantity for a strip geometry in the backgrounds dual to a p-wave and 
to a colorful p + ip superconductor. 

This paper is organized as follows: in section [2] we compute the backreaction of a 3+1 gravity dual to a 
p-wave superconductor in 2+1 dimensions and analyze its thermodynamic properties. In order to compare 
with the colorful superconductor, we review in subsection 12 . 21 the backreaction of the gravity dual of a p + ip 
superconductor. On section [3] we compute the holographic entanglement entropy for a strip geometry in 
both systems, as a function of the temperature and the length of the strip. The conclusions are summarized 
in section |4] 



2 p and p + ip holographic superconductors 

As mentioned in the introduction, the gravity dual to a p-wave superconductor is modeled by an Einstein- 
Yang-Mills (EYM) theory. In [6, 28], the 3+1 dimensional gravity theory dual to a p- wave superconductor 
has been computed in the probe limit. Moreover the authors showed that the p+ip superconductor geometry 
studied in |11| was unstable under small fluctuations, and that the stable configuration was that of the p-wave 
solution. In this section we compute the backreaction of the non-Abelian gauge field on the geometry dual 
to a p-wave superconductor in 3+1 dimensions and compare the results with those for the p + ip case. 

We will work in the simplest set up and consider SU (2) as the gauge group. In the p+ip case, the ansatz for 
the gauge field is such that it breaks the U(l) subgroup of the internal gauge SU(2) and the spatial rotational 
SO(3) group symmetries into a diagonal subgroup of them. Instead, the the p-wave superconductor, breaks 
both U(l) symmetries completely. The gravity solution that describes the strong coupling dynamics of both 
kinds of superconductors is as follows: a charged superconducting layer develops outside the horizon due to 
the interplay between the electric repulsion (with the charged black hole) and the gravitational potential of 
the asymptotically AdS geometry. At high enough temperatures there is no hair outside the black hole and 
the solution is just an AdS-Reissner-Nordstrom (AdSRN) black hole. Below a critical temperature T c a non- 
trivial gauge field with non- vanishing chemical potential on the boundary of the geometry and a sourceless 
non- vanishing condensate in the bulk appears, originating a breaking of the SU(2) gauge symmetry. 
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2.1 p-wave superconductor in 2+1 dimensions 
2.1.1 Solution 



We start from 3 + 1 SU (2) Yang Mills Theory in AdS gravity (see for a review about solutions for this 
theory), the Lagrangian density is 

K f 4} C = R — 2A— ^Tr(F^F^) (3) 

where A = — -j^, is the gravitational constant in four dimensions and the field strength of the SU(2) 
gauge field is written as 

F% = d,Al - d v A; + g YU ^A\Al (4) 

with g YM = ^t^- the parameter that measures the backreaction and g YM the usual Yang-Mills coupling. We 
use latin letters for SU(2) indexes and greek letters for the space-time coordinates. By scaling the gauge 
field as A = we see that the large g YM limit corresponds to the probe (non-backreacting) limit of the 
gauge field. Roughly one can think that ^ — counts the degrees of freedom of the dual field theory that 
are charged under the SU(2) gauge group. Moreover, -j— counts the total number of degrees of freedom. 

K (4) 

Considering backreaction of the gauge filed amounts to say that the number of charged states is of the same 
order as the number of degrees of freedom of the system. 
The equations of motion following from the action are 

G,.v=R^-\g^R = ^g„ + \Tr[F^]^ 9 -fTr[F lp F^P] (5) 



(6) 



we propose the following ansatz [9l [30] 

ds 2 = -M(r)a(r) 2 dt 2 + 



dr 2 + r 2 h{r) 2 dx 2 + r 2 h(r)- 2 dy 2 



for the background geometry, and 



A 



M(r) 
(r)r 3 dt + CLi(r)r 1 dx 



(7) 
(8) 

and a a the usual 



for the gauge field. Here we use the matrix-valued notation A = A®T a dx>* with t° = % 
Pauli matrices, the SU(2) generators satisfy [r a ,T b ] = e abc T c . A solution developing u> ^ in the gauge 
field ansatz © breaks the U(l) gauge symmetry associated with rotations around r 3 (usually called 1/(1)3) 
and a h ^ in the metric breaks U(l) xy symmetry associated to rotations on the xy plane. At high enough 
temperatures we expect no hair outside the black hole and the solution with no condensate is AdSRN with 



tu(r) 


= 0, 




h(r) 


= 1 




a(r) 


= 1, 




<j>(r) 


- m( 




M(r) 


= r 2 


2 2 




r 2 



r 



(9) 
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Replacing the ansatz inyo the EYM equations of motion results into five equations, three of them are 
second order differential equations, and the remaining two are first order constraints 



E? 8cr 2 V rh 2 M r J \r h 2 8r/i 2 

JL ( + + • 

h? \ 8r 8rM 2 h 2 a' 



1 ( _J2 , u>( 2 h ' , M ' 



a 

8r 2 h V ~ ' M 2 a 2 J '" \r ' ~a 

9l M ^ .,[2h' M' a' 
M 2 a 2 \ h Ma 



r 2 h 2 M \r a 



(10) 



This system of equations enjoys four scaling symmetries that become useful when numerically solving it, 
they are 

1. a — > Act, (j) ~^ A</> 

2. lu — > Xuj, h — > Xh 

3. M -> A~ 2 M, a -> Act, g YM -> X~ 1 g YM , R -t XR 

4. M -> X 2 M, r -> Xr, <p -> A0, cj -> Aw 

Using these scaling symmetries we can set i? = r/j = 1 and fix the boundary value of the metric functions 
<t(oo) = h(oo) = 1. The geometry and the gauge field must be regular at the horizon which implies the 
following expansion in the IR (small r) 

M = M l (r - r h ) + M 2 {r - r h ) 2 + . . . 
h = h + h 2 (r - r h ] 2 + . . . 
a = cr + 01 (r - r h ) + cr 2 {r - r h ) 2 + ■■ ■ 
lu = lu + io 2 (r - r h ) 2 + Ld 3 (r - r h ) 3 + ■ ■ . 

cf> = Mr - r h ) + fo(r - r h ) 2 + . . . (11) 
On other hand in the UV (large r) the desired behavior is: 

M = r 2 



2 , Ml , H) 2 + P 2 



r 8r 2 
hi H) 2 , 



w = w,. 

r 



32r 4 



6r 3 



, + 2 + ^ + ... (12) 
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To achieve SSB, we look for solutions where the non-normalizable component vanishes cjq = 0. Standard 
AdS/CFT dictionary instruct us to interpret the boundary and sub-leading values of <f> as the chemical 
potential fi and the charge density p of the dual field theory [3T]. Moreover, the sub-leading coefficient M\ in 
the boundary expansion of gu coincides with the regularized Euclidean on-shell action [3] . The normalizable 
coefficient in ui is dual to the vacuum expectation value of the current (J*) cx lj\ and serves as an order 
parameter for the system. 

Solutions of the system (|10[) depend on the four IR coefficients <f>i, ljq, /in, gq and the backreaction pa- 
rameter g YM . All other coefficients in (fTTj) can be written in terms of those. We proceed to integrate the 
equations of motion numerically out from the horizon using a shooting method in order to get the desired 
asymptotic behavior. We explore the range g YM £ [0.85, 24] and observe that the behavior of the functions 
does not change qualitatively as g YM is varied. In figure Q] and [2] we give the plot of the solutions of (fTt)|) with 
boundary conditions (jf 2[) . We use /i to adimensionalize whenever needed. This means that we are working 
in the grand canonical ensemble. 





Figure 1: The dimensionless metric functions a(r) 
and h{r) for g YM = 2 and T = 0.2312^. 



Figure 2: The dimensionless metric function M(r) 
and the gauge field functions w(r) and 4>{r) for 
g YM =2,T = 0.2312 M . 



2.1.2 Thermodynamics 

In this section we compute the thermodynamic quantities associated with the solutions. As we shall see from 



the study of the potential function in the grand canonical ensemblqfl we have a second order phase transition 
between a superconducting and normal symmetric phases. 

The temperature of the dual theory is given by the Hawking temperature of the black hole 

where the second equality comes from the consistency of the series expansion (jf f p that relates the coefficient 



4 To go from the grand canonical ensemble (fixed fi) with free energy Q, to the canonical ensemble (fixed p) with free energy 
F, we should add a boundary term to the Euclidean action. This changes the variational problem and implies the known Gibbs 
relation F = SI + fip. 
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Mi with do and <p\ . The area of the horizon, Ah , gives the entropy 



2n 



-A h 



2ir 2 VT 2 



12 2 



4) (24<rg-0?) 2 



(14) 



where V = J dxdy. In hgure [3] we plot the order parameter cjj(i.e. the VEV of the current (J^)) as a 
function of the temperature. Note that at T = T c the condensate vanishes showing the disappearance of 
the superconducting state for T > T c . From our numerical results we find (Jf) oc (1 - ^-f' 2 near T c and 
therefore the critical exponent takes the value 1/2. Figure |4] shows the behavior of the Bekenstein-Hawking 
entropy (| 14[) as function of the temperature for our solution and the AdSRN black hole. 




Figure 3: The plot shows the normalizable coefficient of the w function which is proportional to the condensate 
(J^). The black, green and blue lines refers to solutions with g YM = 1, 1.5, 2 and T c = 0.0749, 0.1565, 0.2312 
respectively. Note that the condensate vanishes for T > T c . 





Figure 4: The entropy as a function of the tem- 
perature. The blue line is for the superconducting 
phase with g YM = 2 and the red line for the normal 
phase (AdSRN geometry) . There is a second order 
phase transition at T = T c = 0.2312. 



Figure 5: The potential function computed from 
(|2"2"|) as a function of T for g YM = 2. The red line 
is the potential for the RN solution and the blue 
line is for the superconductor case. 



G 



The gauge/gravity correspondence identifies the Euclidean on-shell gravity action Se times the tempera- 
ture T as the grand canonical potential function D, of the system. To compute it we continue to Euclidean 
signature, time being compactified with period ^ to avoid singularities. The on-shell action has a factor y 
due to time integration, writing S on - s heii = Sb Z lk one has 



Sbu 



Ik 



dx dy dry'—g C 



(15) 



where the lagrangian density is given by The yy component of the stress tensor is proportional to the 
metric and then the Einstein equations implie that 



G 



uu 



r 



R 



Then we have 



and from this we obtain 



rni 



~R — G r ~\~ Gf -4- G T 



2 l K (4) £ 



R 



r 2 an 2 {A) 



r 3 Ma ( h 



(16) 



(17) 



(18) 



where ' denotes derivative with respect to the holographic coordinate r. Then, the bulk contribution to the 
on shell action (TPol can be written as 



5, 



hid k 



dx dy dr^ 



"S£ = --2- 

K (4) 



"■Ma 



(19) 



where is the boundary of the space. As usual, in order to have a well defined variational problem when 
imposing Dirichlet boundary conditions on the metric we need to add to the action a Gibbons-Hawking term 



If V 

Sgh = t - / dx dyV~9ao V M n^ = ^- r 2 a 



"(4) 



' (4) 



M' „,/V 2 

— + M I h - 

2 V a r 



(20) 



where n^dx. 



M dr is the outward pointing unit normal vector to the boundary and ^ is the determinant 



of the induced metric on the boundary. Precisely at r = r M (|20p diverges and therefore must be regularized 
adding the intrinsic boundary counter-term 



S, 



1 f V 

/ dx dy^-goo = — T - 



"(4) 



"(4) 



Ma 



(21) 



Finally the dual thermodynamic potential Q results 

Q = lim S n- S hell 

Too— >00 

= lim (Sbulk + Sqh 



S c 



(22) 



Upon regularizing the action the potential Q results to coincide with the sub-leading value of the g tt 
component of the background metric i.e. i7 = M\ [4]. We have verified our numerical solution computing 
Q in both ways finding an excellent agreement. In figure [S] we plot the potential (|2"2")l as function of the 
temperature. As we mentioned above a second order phase transition develops at T = T c : the grand potential 
and the entropy are continuous but 5* is not differentiable. Below T c the system is in the superconducting 
phase, as we increase the temperature above T c the AdSRN geometry dominates the free energy, this models 
a transition from a superconducting to a normal phase. 
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2.2 p + ip wave superconductors 

Here, we review the results of and compare them with the results of the previous section. We will find 
that at T = T c the system has a second order phase transition and for all ranges of temperatures the grand 
potential of the p-wave solution found in previous section is lower than that of the p + ip, implying that the 
stable phase of the system is the p-wave phase in accordance with the stability analysis [6]. 

2.2.1 Solution 

The background and gauge field ansatz for model a, p + ip-wave solution are 



ds 2 
A 



-M(r)dt 2 + r 2 h{r) 2 (dx 2 + dy 2 ) + 
>(r)T 3 dt + u (r)(r 1 dbs + T 2 dy). 



dr 2 



(23) 
(24) 



One important difference with the p-wave superconductor of the previous section arises in the choice of the 
gauge field ansatz that now breaks the U(l)3 x U(l) xy into a diagonal combination. The p-wave case fully 
breaks the U(l)s x U(l) xy . This allows us to use a metric ansatz that is totationally symmetric in the 
xy-plane. 

The equations of motion obtained for this ansatz are four second order differential equation plus a first 
order constraint arising from the rr component of the Einstein equations 



h" 
M" 



h 

~2 

3 

M 

2.9' > 



i 

M 
r 



M' (j) 1 ' 



B?M rM 8M Ar 2 h 2 



Af 1 
~AI + r 

2 J.2 " 



r 2 h 2 

,2 



>-h 2 M 



= -^r- + ^r 



3_ 



At 

72 



2cj>' 
I - 



AI 
1 



Arh 2 
M'u' 



M' 



h' 2 
h r 



b! 


M' 




' 17 


3 
4- -( 


h> 2 + 



2 2 

8r 2 hM 



+ 



Ar 2 h 2 



M 2r 2 h 2 
3w 2 1 



M 2r 2 h 2 



AI 



h' 



h' 












4) 


+ M' 


+ h 
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Ah 2 M 

The equations have three scaling symmetries that will help us to numerically solve the system. They are 
1. u — > Acj, h — > Xh 

R -> XR, g YM - 
y, <f> X(f>, r -t Xr 

Th = 1 and the value of h(r) at the boundary to h(oo) = 1. The IR behavior of 



2. M -> X- 2 AI, 4 

3. M X 2 M, h - 
and allows us to set R 



9y M 

A 



8 



these equations are those of a charged black hole 



M 
h 

ui 



M x {r - r h ) + M 2 (r - r h f + . . 
h + hi(r - r h ) + h 2 (r - r h ) 2 - 
loq + wi(r - r h ) + u 2 (r ~ r h ) 2 
<t>i{r - r h ) + <f) 2 (r - r h ) 2 + ... 



(25) 



where as before we impose the Maxwell potential (f> to vanish at the horizon in order to have a well defined 
gauge field in the Euclidean continuation. On the UV we demand 



M 



2h\r 



h\ 
r 



1 + 



(h\f 
48r 4 



Mj -8h\Ml + p 2 + 2{uj b 1 )/3 
r 8r 2 



+ 



(26) 



Note that for SSB we do not allow for a non-normalizable piece in w. As before the scaling symmetries ^ 
allows to fix h\ = 1. In figure [5] we plot the behavior of the solutions and figure [7] shows the order parameter 
( J].) >cx lo\ as function of the temperature. For T = T C both condensates vanish and a second order phase 
transition onsets. Note that the values of the condensate for the p + ip case are lower than those in the 
p-wave case. 




1.5 2.0 2.5 3.0 3.5 4.0 



Figure 6: Behavior of the dimensionless functions 
M(r), h(r) 7 w(r) and 4>{r), plotted for g YM = 2, T = 
0.2312^. 



Figure 7: The dual theory VEV ( J^) oc uj\ as a func- 
tion of temperature for the case of the p-wave (blue 
line) and colorful (orange line) superconductors for 
g YM = 2. Its vanishing for T > T c = 0.2312, sug- 
gesting a phase transition between a superconducting 
and a normal state. 
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2.2.2 Thermodynamics 



The temperature associated to the background solution is proportional to the derivative of the gu component 
of the metric evaluated at the horizon. In this case one has 



T 



Ml 
2tt 



(27) 



The Bekenstein-Hawking formula, that relates the entropy with the area of the black hole horizon, for the 
present case reads 

S = ^A h = ^rlhl (28) 

«( 4 ) K (4) 

Figure [5] shows the entropy in the p + ip (orange line), p-wave (blue line) and RN (red line) cases. 
The grand potential f2, is given by the sub- leading coefficient of the metric function g u , as 



K (4) 



(29) 



This potential is plotted in figure [§1 it clearly shows that for any given temperature the p-wave solution 
(blue) is preferred over the p + ip (orange) state. For T > T c the system is in the normal phase (red) and 
the condensate (shown in fig. [7]) vanishes. 




Figure 8: Bekenstein-Hawking entropy for the RN 
(in red), p + ip (orange) and p-wave (blue) solutions 
with g YM = 2 and T c = 0.2312. There is a second 
order phase transition for both superconductors at 
T = T r . 




Figure 9: Grand Canonical potential as a function of 
the temperature. RN solution (red), p + ip (orange) 
and p-wave (blue). For all range of temperatures 
below T c the p-wave solution is preferred over the 
colorful one (g YM = 2 and T c = 0.2312). 



3 Holographic Entanglement Entropy 

An holographic prescription to compute entanglement entropy (EE) on the AdS^+i gravity dual of a CFT^ 
was given in |12) in terms of minimal surfaces. The entanglement entropy AdS prescription involves subdi- 
viding the system into two regions, A and it's complement B, and find the minimal static d — 1 dimensional 
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Figure 10: Diagram of the stripe shape for the region A used to compute the entanglement entropy. 



surface (at constant time) 7,4 such that its boundary coincides with the boundary of the subsystem A (see 
figure [101). 

The entanglement entropy between the two regions is proposed to be the classical area of 7_4, 

^ = -^[ d^a^\ (30) 

where g^ d is the induced metric on the surface and K 2 d+1 j is the gravitational constant in d+1 dimensions. 

In [32] and [33] the EE was computed respectively for backgrounds dual to confining large N gauge theories 
and for several black holes geometries. In this section we perform this computation for a general background 
and apply it to the p and p + ip-wave superconductors found on section [2] Note that the prescription to deal 
with the EE is very similar to that made in [34]. [35] to compute VEV of Wilson loops. In the last case the 
object being computed is the minimal area of a string that explores the AdS space with its endpoints fixed 
to the boundary. The following discussion follows closely that made in [36]. 

We write the d+1 background metric as 

ds d+1 = -g tt (r)dt 2 + g XtXt (r)dx 2 + g rr (r)dr 2 , i=l...d-l, (31) 

where r is the holographic coordinate. The region of interest consists in the straight belt in the direction Xj 
with width L on the X\ direction. The static embedding belt ansatz is x\ = £i(C)j %j = Cj^ r = r (C)j with 
j = 2, . . . , d— 1. A diffeomorphism invariance in remains, depending on the context it will be fixed either 
as x\ = ( (global embedding) or r = £. The entropy ([30]) is 

Sa = -5— / d C\/gx 2 x 2 (r) . ■■gx d - 1 x d ^ 1 (r)Jg rr (r)r' 2 + g XlXl {r)x? , (32) 
K (d+i) J 

where A = JdC2 • • -d-Cd-i and ' denotes x derivatives. Defining g X x( r ) = 9x 2 x 2 ( r ) ■ ■ ■ 9x d ~ix d ^i {r) and the 
functions 

f 2 ( r ) =3xx( r )fc*i( r )> V 2 (r) = 9 X x( r )9rr{r) (33) 
the entanglement entropy is written 

Sa = / dcJv 2 (r)r' 2 + P(r)x' 2 . (34) 

K (d+i) J 
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By minimization of (|34[) we obtain 

4(0 ■ ± /y ;'«> (35, 

/( r ) v 7 / ( r ) - / M 

where r = r*o is the the minimum value in the holographic coordinate reached by the surface. Depending on 
the background under study this could be the horizon radius or the end of the space-time. Inverting this 
relation we can read the length of the belt in the X! direction 

f 00 dr dXl - 2 f°° d A r) /(ro) (36) 

L dr J ro f(r) y/p(r)-p(r Y 

We now fix the remaining diffeomorphism invariance as X\(C) = this choice has the advantage of providing 
a complete parametrization of the embedding r(x\), {x\ <G [ — ■§! "§] an d the boundary condition are J"(±|f) = 
oo). From (|35|l in (j34|l the entanglement entropy reads 

SAM -2™-r+ ,, (37) 



«(d+i) J„, v//(r) 2 - W« 



Expression (f3"T|) diverges at r = oo due to the infinite extension of the surface. The interpretation of this 
divergence is that another solution exists, with the same boundary conditions, consisting on two disconnected 
surfaces expanding all along the radial direction. Its area is 

O A />OC 

S AdiBa - 2- / dr r,(r), (38) 

K (d+1) Jr mm 

here r m i n is the minimum value of r allowed for the geometry The EE is defined therefore with respect to 
the reference state 



In what follows we are going to study the EE for the solutions of (|10p and (|2.2.ip . 
In the p-wave case the relevant functions are: 

r 2 

fp(r) = g yy gxx = r 4 , r)l(r) = g yy g rr = (40) 

where the sub-index p reminds that they correspond to the p-wave superconductor. With this, we can 
compute explicitly the quantity (|39p 

( / dr—=^== r dr-4= ] (41) 




/iv^x/r 4 - r 4 Jr min hVN 



On figure [TT] we plot AS a as a function of the length of the strip L. This shows ASj. for different values 
of the backreaction parameter and different values of the temperature. As expected the bottom line is the 
one which has the lowest temperature because as we lower the temperature we must have more degrees of 
freedom that condense. The linear behavior for large values of [iL is a manifestation of the area law proposed 
in (|50|) . On figure [T3] we plot the EE of the condensed (blue line) and normal (red line) phases as a function 
of the temperature and for a constant value the length of the belt L. 



12 



In order to deal with a finite entropy, avoiding the substraction of the disconnected solution, we can write 
the EE as 



<5^(r ) 



4-kA 



K 2 

"-(4) J ro 



77 



dr- 



47rA 



-R 



(42) 



where Sa has dimensions of length with no divergences. The figure [JJJ] shows that the EE for the superconduc- 
tor (blue line) is lowest that for the RN (red line) solution. This is expected because in the superconducting 
state there are condensed degrees of freedom. 

«?41 as* 




An fi A 




Figure 11: Entanglement Entropy as a func- 
tion of the size of the strip for the p-wave so- 
lution. The black, green and blue lines are for 
values g YM = 1,T = 0.0749p, g YM = 1.5, T = 
0.1565p and g YM = 2, T = 0.2312/1 respec- 
tively. 



Figure 12: Entanglement Entropy as a func- 
tion of the size of the strip for the p + ip so- 
lution. The black, green and blue lines are for 
values g YM = 1,T = 0.0749//, g YM = 1.5, T = 
0.1565/t and g YM = 2,T = 0.2312/t respec- 
tively. 



Performing the same analysis for the solutions of (|2.2.1[) we obtain 



fp+ip( r ) 



9yy9xx 



r 4 h\ 



Vp+ip( r ) — 9yy9rr — 



M 



and the following EE: 



AS A 




dr- 



r 3 h 3 



dr 



rh 



(43) 



(44) 



In figure [12] we show the behavior of AS a in this case, and we can perform the same analysis as for the 
p-wave superconductor. Again, a different approach to obtain a non-divergent entropy, as in the previous 
case, avoiding the substraction of the disconnected surface, consists in separate the divergent piece of the 
integral (|37|) and take in account the finite part of it, Sa- In this case: 



<Sa(t- ) 



47rA 



77 



dr- 



! /i 3 



= S A 



AttA 



-K. 



(45) 



"(4) Jr VM^r^h* - r 4 h(r ) 4 <« (4) 

On figure [TBI we plot this finite part (orange line) for g YM = 2 and T c = 0.2312 and show that, as expected, 
is lowest than the EE for the RN solution. Moreover, the figure shows that the EE in the p-wave case is 
lower than in the p + ip superconductor. This suggest that for a given temperature there are more condensed 
degrees of freedom on a p-wave superconductor. 
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Figure 13: Entanglement Entropy as a function of T and fixed fiL = 3, g YM — 2 and T c — 0.2312. The 
blue line is for the p wave solution, the orange line is for the p + ip superconductor and the red one is for 
the Reissner-Nordstrom solution. 

4 Summary 

On this work we studied the holographic dual to a p-wave and a p + ip-wave superconductors in 3+1 di- 
mensions. We computed the backreaction of the p-wave solution and studied its thermodynamics properties. 
As expected, and in contrast with the solution in 4+1 dimensions studied in [9], we found a second order 
phase transition between the normal and superconducting phases. Later, we reviewed the backreaction of 
the gauge field on the geometry of the colorful black hole and compared it with our solution. From the study 
of the thermodynamic quantities and in particular from its grand canonical potential, we noted that for a 
fixed value of the temperature the p-wave solution has less potential and then it is preferred. We related 
this with the fact that the p + ip solution is unstable under small fluctuations and it decays into the p-wave 
background. 

Finally, using the holographic proposal given in |12| to compute the entanglement entropy of a quantum 
field theory using its gravity dual, we computed it for both solutions on a straight belt geometry as a function 
of the temperature and of the size of the belt. We observed that for both cases the EE behaves linearly for 
large values of L, which confirms the proposed area law. The EE vs L plots moves to large values of ASj\ 
as we increase the temperature. As a function of the temperature we observe that the largest EE is for the 
RN solution. This is expected because the superconductor has condensed degrees of freedom. Moreover the 
p-wave solution presents more condensed degrees of freedom than the colorful black hole. That may explain 
the fact that the condensate value is larger for the p-wave system at a given temperature. 
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